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Abstract 

Large N geometric transitions and the Dijkgraaf-Vafa conjecture 
suggest a deep relationship between the sum over planar diagrams 
and Calabi-Yau threefolds. We explore this correspondence in details, 
explaining how to construct the Calabi-Yau for a large class of M- 
matrix models, and how the geometry encodes the correlators. We 
engineer in particular two-matrix theories with potentials Y) 
that reduce to arbitrary functions in the commutative limit. We apply 
the method to calculate all correlators {ti XP) and (trYP) in mod- 
els of the form W{X,Y) = V{X) + U{Y) - XY and W{X,Y) ^ 
V{X) + YU{Y^) + XY^. The solution of the latter example was 
not known, but when [/ is a constant we are able to solve the loop 
equations, finding a precise match with the geometric approach. We 
also discuss special geometry in multi-matrix models, and we derive an 
important property, the entanglement of eigenvalues, governing the ex- 
pansion around classical vacua for which the matrices do not commute. 
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1 Introduction 

Ever since the work of 't Hooft on large n QCD [1], finding methods to sum 
up planar diagrams has remained a central theme in mathematical physics. 
Yet, even in the best studied and simplest case of zero-dimensional integrals 
over hermitian n x n matrices, only very few tcchniqTics arc available [2]. 
The usual approach is to solve a set of saddle point equations by making 
an ansatz for the analytic structure of the solution [3]. The saddle point 
equations are themselves derived after reducing the number of degrees of 
freedom from ~ to ~ ra by integrating over angular variables, which can 
necessitate sophisticated methods [4, 5, 6] that apply only in special cases. 
Another approach is to write down Schwinger-Dyson equations, called loop 
equations in this context [7]. This method is very general and does not use 
assumptions on the analytic structure of the solution. However, it is usually 
very difficult to find a finite set of equations that closes under the correlators 
one wish to calculate. This scarcity of tools and of solvable examples makes 
the search for alternative strategies particularly useful. 

It is well-known that the solution of the one-matrix model can be ex- 
pressed in terms of a complex algebraic curve (see for example [8] or the 
appendix of [9] for reviews). The idea we wish to pursue in the following is 
that for multi-matrix models it may be useful to consider higher dimensional 
Calabi-Yau spaces. This is suggested by conjectures by Vafa and collabora- 
tors on the strongly coupled dynamics of four dimensional Af = 1 supersym- 
metric gauge theories [10, 11, 12, 13]. We will engineer multi-matrix models 
with potentials of the form 



where E{z, w) = J2t^oo -^i(^) ^* be expanded in terms of entire func- 
tions Ei and Co is a small contour encircling z = 0. This is a large and 
interesting class of models. In the one- or two-matrix cases, the commuta- 
tive limit of (1.1) can be an arbitrary function. 

The main body of the paper is divided into four parts. In Section 2 
and Appendix A, we discuss special geometry relations for multi-matrix in- 
tegrals. Special geometry is automatically built in the Calabi-Yau geome- 
try, and it is important to have a satisfying understanding from the matrix 
model perspective as well. Deriving the relations, we uncover new inter- 
esting properties of the saddle point equations and their solutions (eigen- 
value entanglement) around vacua where the matrices do not commute. 
These developments in standard matrix model technology are not strictly 




(1.1) 
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required to understand the geometric approach, that we describe in de- 
tails in Sections 3 and 4. The method is apphed in particular to com- 
pute the resolvents of both matrices X and Y for the standard [14] two- 
matrix model W{X, Y) = V{X) + U{Y) - XY as well as for the model 
W{X,Y) = + v{X) + YUiY"^), for arbitrary polynomials V and U. 

In Section 5, we use the loop equation technique to solve for all correla- 
tors {ti X^Y'i) in the latter example when J7 is a constant, finding a perfect 
match with the geometric approach. Finally, we discuss open directions of 
research in Section 6. 



2 The geometry of matrix models 

The existence of special geometry relations constraining the sum over planar 
diagrams is a basic ingredient of the geometric approach to matrix models. 
This is an essentially unexplored subject in the case of multi-matrix integrals, 
and thus we provide a rather detailed, albeit incomplete, discussion in the 
present Section and in Appendix A. 

A matrix model with M matrices has several one-matrix model effective 
descriptions obtained by integrating out all but one matrix, for example 



By analogy with the ordinary one-matrix model with polynomial potential, 

it is natural to suspect that special geometry relations could be derived for 
the multi-matrix models by working with the effective descriptions. The 
difficulty is that the effective potential V^s can have a non-trivial analytic 
structure. The main subtleties come from the classical solutions for which 
the matrices do not commute. By studying the properties of the effective 
potentials, we derive a remarkable property of the eigenvalue distributions 
governing the expansion around such solutions. Full details are provided in 
Appendix A for two models that display some generic features. 

The various sets of special geometry relations, derived for each well- 
behaved one-matrix effective description, must be equivalent. In the ge- 
ometric approach, this non-trivial constraint is automatically satisfied by 
expressing the solution in terms of integrals over a higher dimensional com- 
plex variety equipped with prefered coordinates associated with the various 
matrices. 




i=2 



(2.1) 
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2.1 Definitions 

We consider a model with M hermitian n x n matrices Xj- and polynomial 
potential W{Xi, . . . Xm)- The classical equations of motion 

A: tidW{Xi,...,XM) = (2.2) 

define an algebra A with a set of D/<:-dimensional irreducible representations 
Rk- In the generic case there are no moduli and the index K is discrete. 

A particular solution to (2.2) of the form 

R = ^ukRk , ^ riKDK = n , (2.3) 

K 

is called a classical vacuum. The partition function J^{R) is defined by the 
formula 

^n^myS' = r dXi ■ ■ ■ dXM W^(Xi,...,Xm) ^ (2.4) 

•J R— planar 

where the subscript "i?-planar" in (2.4) means that we sum up only the 
contributions from the planar diagrams in the perturbative Feynman expan- 
sion around the classical solution R. The partition function depends of the 
integers uk in (2.3), or more conveniently on the variables 

Sk = ^S (2.5) 
n 

that satisfy the identity 

Y,DkSk = S. (2.6) 

K 

In general, the perturbative expansion defining involves ghost matrices 
as explained in [15]. The coupling S plays the role of h. An important 
result is that the sum over planar diagrams is always absolutely convergent 
for small enough S [16]. For that reason the convergence properties of the 
integral J dXi ■ ■ ■ dXM exp(— n tr W/ S) arc irrelevant. It can always be made 
convergent without loss of generality by performing an analytic continuation 
and/or turning on some couplings. If the integral is convergent, then its strict 
large n limit yields T{Rram) for the representation .Rmin that minimizes W 
(assuming that it is unique). A more useful statement is that for general 
R, J'{R) is the large n limit of the formal perturbative series around the 
classical solution R. In this limit the variables Sk defined in (2.5) become 
continuous. 
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Correlators are defined by 

{0)r = e-"'^(«)/^' [ dXi • • • d^M Oe"t ^^-^(^i'-'^-) . (2.7) 

J ii— planar 

We will not systematically indicate the dependence in R. Basic objects that 
we want to compute are the resolvents of the matrices Xjn defined on the 
complex plane by 

The resolvents are the generating functions of the correlators {tr Xm) for 
arbitrary positive integers p. They are analytic functions with square root 
branch cuts. The physical sheet is defined by the asymptotics 

The densities of eigenvalues p^"" of the matrices are read off from dis- 
continuities of the resolvents across branch cuts on the physical sheet, 



27rS 

and we have 

9 



P^^ (a;) = 7^ {g^^ {x + ie) - g^^ {x - ie)) , (2.10) 



r-PM£)d.. ^^^^^ 

J ~oc X Z 



2.2 Special geometry 
2.2.1 Generalities 

We work with the one-matrix model whose potential is defined by (2.1). For 
a general vacuum (2.3), the eigenvalues of X form cuts I§ I < k < Dk, 
that generically do not overlap. The number of eigenvalues in the cuts, 
or equivalently the variables Sk defined in (2.5), are expressed in terms of 
period integrals, 

5k = ^ / g''{x)dx, (2.12) 

where the contours ^. encircle the cuts ^ on the physical sheet for g-^ 
(see Figure 1 for the definition of various contours). The formula (2.12) 
could be written in many different-looking but equivalent ways, using other 
matrices Xj of the original multi-matrix model. Taking into account the 



624 



PLANAR DIAGRAMS AND CALABI-YAU SPACES 




Figure 1: Definition of the contours a^, I^k K' used in tlie main 

text (similar contours like or are also used). The eigenvalues that 
classically sit in the representations Rk and Rk' are spread along the cuts 
~ ['^Kibx] and I^, = [ax'-ibx'] respectively. 



constraints (2.12), the eigenvalues adjust to make the partition function 

T = -5V;ff(xi, . . . ,x„) + — ^ln|a;i - Xj\ (2.13) 



extremal. The corresponding saddle point equations 

E:r^ = o (2.14) 



5Feff , 25 

-n 



are interpreted as vanishing force conditions. The total force on a test eigen- 
value on the complex plane 

f{x) = Mx) + 2g''{x) (2.15) 

is the sum of a background force 

Mx) = -n^ , (2.16) 

Xi=X 

and of the Coulomb term 2g^[x). Since varying Sk amounts to moving 
eigenvalues from one cut to another [13], we expect that dT/dSx will be 
expressed in terms of countour integrals of the differential fdx. If the term 
/bdx docs not contribute to the integrals, as in the ordinary one-matrix 
model for which it is an exact differential — diy, the special geometry rela- 
tions associated with (2.12) straightforwardly follow, for example 



(p 5^ (x) da; -I- counterterm , (2-17) 



dSK 

where the non-compact 7-cycles are defined in Figure 1. However, in general, 
/bdx do contribute, and a more detailed discussion is needed. 
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2.2.2 The commutative case 

There is a class of models for which the classical equations of motion imply 
that the matrices commute with each other. The non-commutative structure 
can then, in some sense, be neglected. This is also true more generally for 
vacua built from one dimensional representations only. A typical example is 
the much-studied two-matrix model with potential 

W{X,Y) = V{X) + U{Y)-XY, (2.18) 

where V and U arc polynomials of degree dx + 1 and -|- 1 respectively. 
This model was used for example in the study of two dimensional gravity 
coupled to c < 1 matter [8]. The classical equations of motion 

A: X = U'{Y) = Q{Y) , Y = V'{X) = P{X) , (2.19) 

imply that [X, Y] = 0, and thus the only irreducible representations are the 
dxdy one dimensional representations given hy X = xl and Y = yl with 

x = QoP{x), y = P{x) ^ y = PoQ{y), x = Qiy). (2.20) 

The effective potential for X is the sum of a classical part (tr V) jn and of 
a quantum term proportional to the partition function of the one-matrix 
model in an external field. For cubic this partition function is related 
to the Kontsevich integral for two-dimensional topological gravity [17], and 
was studied in [18]. In general, a qualitative picture of the analytic structure 
of the background force can be obtained by assuming that S is very small. 
In this limit the matrix Y can be integrated out classically and we get 

/b,ci(.T) = -P(x) + Q-nx). (2.21) 

The fTinctional inverse Q^^ of the polynomial Q is a multi- valued function 
with dy sheets. There are — 1 branching points at a; = Q'~^(0). The 
analytic structure of /b at finite 5 is a perturbation of the analytic structure 
of /b,cl) with the same number of sheets and of branching points. This can 
be explicitly checked when \J is cubic by using the results of [18], and no 
other singularities are expected. In particular, does not have branch cuts 
on the support of p^, and thus does not contribute to contour integrals over 
a-, /9- or 7-cyclcs. The relations (2.17) are thus valid. Moreover, the saddle 
point equation (2.14) has a simple n ^ oo limit, 

h{x) + g^{x + ie)^g^{x-it) = ^ for x G Support[/9^] , (2.22) 

showing that is a [dy + l)-shcctcd analytic function (the physical sheet 
and the dy additional sheets of /b). This is consistent with the fact that 
is an algebraic function of degree dy-l-l (see [14] and Section 4). A symmetric 
discussion applies to the effective description in terms of the matrix Y . 
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2.2.3 The non-commutative case 

Things are much more subtle and interesting when higher dimensional rep- 
resentations are present. The naive classical limit of a saddle point equation 
like (2.22) would yield the one-dimensional representations only. The emer- 
gence of the higher dimensional representations, which is a classical effect in 
the description in terms of the M matrices Xi, must come out as a quan- 
tum effect in the effective description in terms of the single matrix X. This 
means that some of the quantum corrections in /b are singular and can have 
a finite S* ^ limit that modifies the naive classical limit. A related crucial 
fact is that, even though the equations (2.14) are expressed in terms of the 
eigenvalues of X only, they must know about correlations between X and 
the other matrices that do not commute with X. This implies that in the 
continuum limit, (2.14) cannot be expressed in terms of only. Wc will 
argue that the density is still unambiguously determined, because there 
is a constraint, the entanglement of eigenvalues, that relates p^ on the var- 
ious intervals I^j, associated with a given representation Rk- Intuitively, 
this constraint comes from the U(n) gauge symmetry of the original matrix 
model. 

Those facts are very general but are best illustrated on an example. Let 
us consider the model with potential 

W{X, Y) = XY^ + aY + V{X) (2.23) 

for an arbitrary polynomial V of degree d+ 1. It is useful to separate V'{X) 
into an even and an odd part, 

-V'{X) = Fi{X^) + XF2iX^) . (2.24) 

The classical equations of motion 

A: {X,Y} = -a, Y^ = Fi{X^) + XF2iX^), (2.25) 

admit d + 2 irreducible representations Rk of dimension one, X = xl and 
Y = yl with 

2xy = -a, y'' = -V'{x), (2.26) 

and [(d — l)/2] irreducible representations Rk of dimension two for which, 
in terms of Pauli matrices, X = xa^ and Y = yia^ + y^a^ with 

2xy^ = -a, y'' = yl + yl = Fi{x''), F2(x2) = 0. (2.27) 

If only the one dimensional representations are considered, then the saddle 
point equations for the matrix X have standard features and can be solved 
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-xg -X7 -xg -X5 -X4 5 -X2 -xi 

I • III 1 

I x=0 III 1,^ , 

Xi X2 X3 X4 X5 Xg 8+ X7 Xg 

Figure 2: Eigenvalues filling a two dimensional representation Rk for = 4 
(thick lines) and their images with respect to x = (thin lines). When 
the eigenvalues feel a quantum even force that compensates for the 
non-zero classical force Fi — a^/(4x^) in (2.27). 



by the usual tricks. This was done for a = and with only one cut that 

do not intersect the origin in [19]. We are interested in the general case. 
The effective potential for X, the background force and the saddle point 
equations can be calculated exactly, 

= -V'(.) -^-^^-1- = -V'i.) - ^ + , (2.29) 

3 ■' 



-n^,)-^--E^— + -E^— = 0- (2-30) 

4x,- n ^-^ Xi + Xj n ^ Xi — Xj 

It is important to assume that a 7^ 0, because otherwise the vacua X = 0, 
Y = iy/ -Fi(O) a^, for which the Z2 symmetry Y — > —Y is broken, can no 
longer be described after the exact integration over Y. If needed, the limit 
a — > can be taken on the final formulas. The singularities in (2.28) and 
(2.29) can be understood as follows. Classically, integrating out Y amounts 
to solving the equation 

{X,Y} + a = (2.31) 

for a given matrix X, and this implies that [X^,y] = 0. Generically (2.31) 
has a unique solution Y = —aX~^/2, which is singular if X is not invertible, 
yielding poles at Xj = 0. More interestingly, if two eigenvalues of X are 
such that Xi = —xj, then a different class of solutions to (2.31) is possible, 
corresponding to the two-dimensional representations. In that case, the S — > 
limit of the term l/(a:j + Xj) in (2.30) is non-zero and compensates 

for the non-zero even background force in (2.27). 

A typical configuration with a two-dimensional representation Rk, corre- 
sponding to a solution = x'j^ of (2.27), is depicted in Figure 2. We put hk 
eigenvalues xi, . . . , x„^ at x = —X(^k) and nx eigenvalues . . . , X2nK 

at X = X(^x)- Turning on S, the eigenvalues fill intervals = [aKjix] and 
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~ ['^'k^ ^'k\ ^-round x^^) ^^nd — respectively. The attractive force be- 
tween eigenvalues and their images, described by the term 5 \/{xi + xj), 
induces strong correlations between the eigenvalues around +xk and the 
eigenvalues around —xk- If we choose the labels such that xi < X2 < ■ ■ ■ < 
X2ni^, classically xi = • • • = = -Xf^jq and = • • • = X2nx = ^(K)- 

Quantum mechanically, the only stable, finite energy configurations are then 
such that —X2nK < < ~X2nK-i < for 1 < z < riK — 1, or 

Xi < —X2nK-i+i < for 1 < i < uk — 1 and Xn^ < —Xn^+i- This en- 
tanglement of eigenvalues can be proven as follows. If more than one image 
eigenvalue —xj, j > nx, were in a given interval ],Tj,,Tj+i[, i < uk — 1, then, 
because of the attractive force between eigenvalues and image eigenvalues, 
the system would collapse to a singular configuration with divergent en- 
ergy. If there were some intervals without image eigenvalues then, 
due to the repulsive force, the eigenvalues that are not entangled would not 
converge to the equilibrium positions xk or —xk for the two-dimensional 
representations Rk when 5 — > 0. 

A direct consequence of the entanglement is that, in the continuum large 
n limit, the intervals Ik and are symmetric, h'j^ = —clk and = —bx, 
and 

p^{x) = p^(-x) for xei^. (2.32) 

This shows that the background force (2.29) has branch cuts that coincide 
with those of g-^ , and in particular will contribute to contour integrals. 
Moreover, the continuum limit of the saddle point equation (2.14, 2.30) 
must be studied with care. One might want to replace f\){x) by {f-b{x + ie) + 
fh{x — ie))/2, but this is not correct. As explained in Appendix A, Section 
1, if d^{x) and S~{x) are the even functions defined in Figure 2, then the 
continuum limit is 

^ {fh{x + ie) + fb{x - ie)) + {x + ie) + g^{x- ie) = 

5~ ix) 

Sp^{-x) In — f4 ' for X G Supportfp^l . (2.33) 
o~^{x) 

The even "quantum" force on the right hand side of (2.33) can compen- 
sate for the even classical force Fi(.t^) — a^/(4x^) in the S" — > limit. 
Consequently, there can exist classical equilibrium configurations for which 
only the odd part xF2{x'^) of the classical force vanishes. Those configu- 
rations correspond precisely to the two dimensional representations (2.27). 
By analysing (2.32) and (2.33), we further show in Appendix A that g^ is 
generically three-sheeted. It satisfies a degree three algebraic equation that 
we compute by two different methods in Sections 4 and 5. 
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2.2.4 Lesson 

For higher dimensional representations, we have seen that the continuum 
saddle point equation involves additional unknown functions and is supple- 
mented by conditions coming from the entanglement of eigenvalues. Once 
this is understood, one can proceed and try to derive the special geometry 
relations. This is done is Appendix A for the model (2.23), both in the 
description in terms of X discussed above and in the description in terms of 
Y. The same analysis applies straightforwardly in a large class of models, 
including the famous Af = 1* theory, in which case an ansatz corresponding 
to eigenvalue entanglement was first made in [21]. 

In general, the effective potential, as a function of the eigenvalues of a 
given matrix X, has singularities when the eigenvalues satisfy some condi- 
tions corresponding to the emergence of representations for which X does 
not commute with all the other matrices. In our case the condition was 
but other type of conditions, for example Xi — Xj G Z, are possi- 
ble. The logarithmic singularities that we found in (2.28) are singled out by 
the fact that V^s is a well-defined functional of in the continuum limit 
(which is necessary for the idea of a master field, and the corresponding 
factorization of correlation functions, to apply), unlike its derivatives, that 
know about the correlations between X and the other matrices. The gen- 
erated force must be attractive in order to stabilize new higher dimensional 
solutions in the classical limit. These features presumably imply eigenvalue 
entanglement and special geometry in a large class of models. 

2.3 Solving the constraints from special geometry 

A natural way to implement special geometry is to express the solution in 
terms of a Calabi-Yau manifold with nowhere vanishing holomorphic top- 
form Q, 



The Ak and Ck are holomorphic spheres in one-to-one correspondence with 
the (Xk"^ ^^'^ cycles. A very important additional property is that, at 
least in a large class of examples, there are privileged coordinates associated 
with the matrices, or more precisely with the generators of the center of 
the algebra A. Integrating over the coordinates in different orders in (2.34), 
we can obtain the description in terms of the various matrices. We will 
argue in Section 4 that this property makes possible the computation of the 
resolvents. 




(2.34) 
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3 Matrix models and Calabi-Yau threefolds 

Let us consider type IIB string theory one some non-compact Calabi-Yau 
threefold. Let us assume that there is a with normal bundle M and 
quantum volume V in the geometry. By wrapping D5 branes on the 
P^, we engineer a four dimensional A/" = 1 supersymmetric gauge theory, 
with bare Yang-Mills coupling constant 5^^^ ~ 1/^ and gauge group U(A^). 
The deformation space of the P^, or equivalently the space of fluctuations 
of the brane, is generated locally by the holomorphic sections of the normal 
bundle M. Physically, this means that if M has h linearly independent 
holomorphic sections, the gauge theory is coupled to h scalar superfields in 
the adjoint representation. For a generic geometry, there is an obstruction to 
the deformation space, which in favorable cases can be interpreted physically 
in terms of a superpotential W. 

If /i < 2, the resulting gauge theory is asymptotically free and dynami- 
cally generates a mass gap A. The renormalized Yang-Mills coupling gj^ at 
scale ji grows in the IR, and becomes infinite at some scale /Xc ~ A. Beyond 
that scale the fundamental gauge fields strongly fluctuate and are no longer 
appropriate degrees of freedom. Prom the type IIB perspective the RG flow 
corresponds to the shrinking of the P^ to zero size. At scale n = /ic, the P-*^ 
vanishes and the original smooth Calabi-Yau is replaced by a singular geom- 
etry. What happens for /i < Hc has not been derived from first principles, 
but there are remarkable conjectures that partially address the problem. On 
the gauge theory side, it is believed that good degrees of freedom are given 
by the so-called glucball superfields Si, and that the exact quantum super- 
potential VFgiuebaiis ^oi the Si can be calculated from the partition function of 
an associated h-matrix model whose potential is equal to the superpotential 
of the gauge theory [13]. On the string theory side, it is believed that the 
singular Calabi-Yau is deformed to a smooth Calabi-Yau. The branes and 
the P^ disappear and are replaced by three-form fiux through three-spheres. 
The non-zero flux generates a superpotential Wflux that can be computed 
from special geometry [22]. The consistency between the conjectures, 

W^glueballs = W^flux , (3.1) 

implies a highly non-trivial relationship between the partition function of the 
matrix model and the Calabi-Yau geometry, whose form coincide precisely 
with (2.34). This relationship is indicated by a question mark in Figure 3. 

The great advantage of this gauge theory/string theory set-up is that we 
have a precise recipe, explained below, to construct the relevant Calabi-Yau 
spaces. On the other hand, it is important to emphasize that there is no 
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Geometric transition conjecture 



Type IIB D5 branes 
wrapped on PICCY3, 
size(Pl) = V 



Blow down 



Singular CY3 
V=0 



Deformation 



Three-form flux 
through S^CCYa 
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N=l super Yang-Mills 
with adjoint scalars, 
coupling l/g2 ~ V 



RGflow 



Strong coupling 




RGflow 



Glueballs S with 

W= VsF 
F = matrix model 
partition function 



Dijkgraaf-Vafa conjecture 



Figure 3: Chain of conjectured dualities suggesting a relation between 
Calabi-Yau spaces and matrix models. The question mark emphasizes the 
link that we want to explore. See details in the main text. 

general proof of the validity of the approach. For example, the full relation- 
ship between J\f = 1 gauge theories and matrix models can be proven only 
in special cases and under some assumptions [23]. A perturbative argument 
can also be given [24], but it does not apply to the objects that are relevant 
for us. In particular, there is no understanding of A/" = 1 special geometry 
[25] in gauge theory from first principles. 

As is clear from Figure 3, there are three geometries that play a role, and 
that we discuss in turn in the following. The first geometry is the original 
Calabi-Yau space M. containing the P^s with normal bunble J\f on which 
we wrap the D5 brancs. This space A4 is called the resolved Calabi-Yau, 
because it is a small resolution of the singular Calabi-Yau A4o obtained from 
A4 by shrinking the (the resolution is small because the singular points are 
replaced by curves). Precisely, there exists a blow down map tt : M. ^ Mq 
that locally induces a birational isomorphism between A^\{P^} and Aio\{p} 
where p is a singularity and ir^^ip) = P^. Finally there is the smooth 
deformed space Ai obtained by perturbing the algebraic equation defining 
A^o- All we need to know about M is this algebraic equation, that encodes 
the complex structure, because the metric datas are irrelevant in integrals 



like (2.34). 
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3.1 The resolved geometry 
3.1.1 The transition functions 

Let us describe M. by two coordinate patches {z^wi,W2) and [z' ^w'-^^w'^), 
where z and z' = 1/ z are the stereographic coordinates on the two-spheres, 
and consider the following transition functions 

z' = l/z, w[ = z-^'wi , W2 = z-'^Wi . (3.2) 

If n > and m < (other cases could be discussed as well, but arc irrelevant 
for our purposes), the geometry (3.2) has a (n + l)-dimensional continuous 
family of P-'^s that sit at 

n+1 

wi{z) = Xiz'-^ , W2{z) = . (3.3) 
1=1 

The normal bundle to the P-'^s is M = 0{n) © 0{m) by construction. The 
functions wi{z) and W2{z) in (3.3) are characterized by the fact that they 
must define globally holomorphic sections of J\f. This is indeed the case 
because in the other coordinate patch, (3.2) yields 

n+1 

■^'i (^') = X] ' W2{z') = . (3.4) 

i=l 

The parameters Xi, 1 < i < n + 1, span the versal deformation space of the 
P^s. A more interesting geometry is obtained by perturbing the transition 
functions given by (3.2). The perturbation is described in terms of a geomet- 
ric potential. The potential is a function E(z, w) of two complex variables 
that can be Laurent expanded in terms of entire functions Ei, 

+ 00 

E{z,w) = Ei{w)z\ (3.5) 

i=—oo 

The most general geometry we will consider is given by 

z' = l/z, w[ = z-""wi, w'2 = z-"'w2 + dyjE{z,wi). (3.6) 

Since the relation between wi and w[ is unchanged with respect to (3.2), 
the most general holomorphic section {wi{z) , W2{z)) of J\f is still such that 

n+1 n+1 

^il-^) = Xiz""'^ <S=4> w'i{z') = X ajjz'""*"*"^ . (3.7) 
i=i i=i 
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On the other hand, W2{z) must now be adjusted non-trivially to insure that 
W2{z') is holomorphic. Because of the term z~™"W2 = z' '™'u;2 in (3.6), one 
can always choose w^iz) to cancel all the poles in z' ^ for j > |m|, but there 
remains |m| — 1 singular terms that cancel only if the parameters Xi in (3.7) 
satisfy |m| — 1 constraints. Thus we find that the vcrsal deformation space 
of is spanned by n + 1 parameters satisfying |m| — 1 constraints (the same 
result could be straightforwardly obtained for the most general perturbation 
of (3.2), see [26]). Since the tangent bundle to is 0(2), the Calabi-Yau 
condition of vanishing first Chern class for the total space is 



and is thus equivalent to the equality of the number of parameters and the 
number of constraints. 

3.1.2 The superpotential 

The special form (3.6) of the perturbation is particularly interesting because 
in this case the constraints arc integrable and are equivalent to the extrem- 
ization dVF = of a superpotential W{x\, . . . ,x„+i). This was proven in 
[27] in the case where £^(z, w) is regular at 2; = 0, and the same argument 
could be straightforwardly generalized to (3.5). We can actually derive an 
explicit formula for W that makes this property manifest. By using (3.6), 
(3.7) and (3.8) we obtain 



Let us define the contour Cz to encircle both the points = and u = z. If 
we choose in a manifestly z-holomorphic way to be 



M = n + 1 = -m - 1 



(3.8) 



M 




(3.9) 




(3.10) 



and impose the M integrable constraints 




then we 



;et a manifestly z'-holomorphic w'^ 




(3.12) 
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The superpotential is 

W{xi,...,xm) = ^ z-^-^E(^z.Y,j.,z'-^)dz. (3.13) 

For example, in the case of one variable we obtain W{x) = Ei{x), and in 
the case of two variables 

W{x,y) = Y,E^liW/i[. (3.14) 

i>0 

It might seem that different geometries (3.6) can yield the same W, but this 
is not the case. For example, in the case M = 2, the holomorphic change of 
variables 

W2 ^ W2 - ^ E2_i{wi)z~'~'^ , 
i<0 

i-2 (3.15) 

i>2 j=0 

put the geometry in the form 

z' = 1/z , w[ = Z~^Wl , W2 = Z^W2 + ^ E'2_i{wi)z'^~' , (3.16) 

i>0 

with the conditions 

^({_).(0) = for 0<j<i-l (3.17) 

on the derivatives of the functions -B^. There is thus a unique geometry for 
a given superpotential. 

3.1.3 The non-commutative structure 

Equation (3.14) shows that we can engineer an arbitrary superpotential of 
two variables when only one branc is wrapped on P^. However, ordering 
ambiguities arise when N > 2 D5 branes are considered. Eventually we 
expect that only special, in some sense integrable, matrix models can be 
constructed. The rigorous computation of the matrix model potential W as 
a function of matrices Xi would require an extensive use of the theory of 
branes in Calabi-Yau threcfolds, and this is beyong the scope of the present 
paper. There is however little doubt as to what the correct answer must be. 
Our proposal is that the parameters Xi parametrizing the fluctuations of the 
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brancs in (3.7) must be replaced by hermitian matrices Xj. Supersymmetry 
then implies (3.11) with xj Xj. Remarkably, with this prescription for 
the ordering, the constraints are still integrable, and the potential is given 
by (1.1). 

In the two-matrix case, on which we shall mainly focus in the following, 
the above prescription tells that the non-commutative version of a monomial 
x^y^ is simply obtained by expanding tvlX + Yy~^^ and picking up the terms 
with the right degrees in X and Y. For example, x'^y^ — > tr(2X^y^/3 + 
XYXY/3). The general formula is 



xY — ' , \, z-^-^ \x(X + yz)'+^ . (3.18) 

For polynomials Vx, Vy and F, the geometric potential 

E{z,w) = z^Vx{w) + z^Vy{w/z) + ^^jjl (3.19) 

yields models of the form 

W{X, Y) = VxiX) + VriY) + V{X + Y) . (3.20) 

The term (1 — l/z)~^ in (3.19) must be understood as a formal power series 
in z~^, for which only the first dcgT^ + 1 terms are relevant. Special cases 
of (3.20) have been obtained by considering Calabi-Yau threefolds that are 
monodromic fibrations of ALE spaces with ADE singularities [12]. The A- 
series correspond to degV = 2, the D-series to degF = degT^ = 3 and the 
E-series to degV = 3, degVy = 4 and 4 < degVx < 6. For those models, 
the deformed geometry M. can in principle be calculated by using the results 
of [28]. 



3.2 Blowing down 
3.2.1 Theorems 

Let us first state an important theorem by Laufer [29]: 

Theorem (Laufer): Let Mq be an analytic space of dimension D > 3 with 
an isolated singularity at p. Suppose that there exists a non-zero holomor- 
phic D-form 0, on A^o\{p}- Let tt : A4 ^ be a resolution of Aio- Sup- 
pose that the exceptional set A = Tr~^{p) is one-dimensional and irreducible. 
Then A is isomorphic to and D = 3. Moreover, the normal bundle of 
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in M must be either N = 0{-l) © C(-l), or A/" = 0(0) © 0(-2), or 

It is known in general that A/" is a direct sum of hne bundles [30], M = 
0{n) © 0{m). The condition on the normal bundle in Theorem 1 is thus 
equivalent to asymptotic freedom, M < 2. This is consistent with the fact 
that asymptotically free theories are the one for which we can expect the 
to be exceptional. The situation for M > 3 is more subtle. An important 
point is that the normal bundle to the in the geometry (3.6) is changed 
when the perturbation E is added. From the gauge theory point of view, 
turning on the superpotential amounts to giving a mass to the chiral multi- 
plets, and gencrically the theory flows to the pure gauge theory in the IR. 
More precisely, the number of massless chiral multiplcts in a given vacuum is 
equal to the corank of the Hessian of W at the corresponding critical point. 
This translates mathematically in the following 

Conjecture (RG flow): Consider the geometry (3.6) for m = — n — 2 and 
associated superpotential W given by (3.13). Let M be the normal bundle 
of a P^ that sits at a given critical point of W . Let r be the corank of the 
Hessian of W at the critical point. Then N = 0{r - 1) ® 0{-r - 1). 

We give an elementary proof of this conjecture for n = 1 in the Appendix. 
This result suggests that models with an arbitrary number of multiplets, or 
integrals over an arbitrary number of matrices, may still be described by a 
geometric transition as in Figure 3. The difhculty is that Laufer's theorem 
implies that the blow down map must be singular when the parameters are 
adjusted in such a way that the corank of the Hessian jumps to a value 
greater than two. Because of this complication, we restrict ourselves to 
two-matrix models in the following. 



3.2.2 Examples 

The blow down map : M. Mq is given explicitly by four globally 
holomorphic functions tTj, 1 < z < 4. The mapping tt must be a birational 
isomorphism except on the P^s that are mapped onto the singular points of 
A^o- We present the construction for three examples. 

Example 1: The resolved geometry 

z' = 1/ z , w'l = wi , W2 = z^W2 + zP{wi) , (3-21) 
engineers the one-matrix model with a potential W{X) such that W' = P. 
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The P^s sit at 

wi {z)=x = w[ {z') , W2{z) = = w'2 {z') , (3.22) 

with P(x) = 0. The blow down is straightforwardly found to be 

TTi = u;i = w[ , (3.23) 

7r2 = 2zw2 + Piwi) = 2z'w'2 - Piw[) , (3.24) 

TTg = - zP{wi) - z^W2 = -w'2 - z' P{w\) + z'^^a , (3.25) 

7r4 = u;2 + zP(i(;i) + z^u;2 = - '^'^'('i^i) + -^'^"^2 ■ (3-26) 

As required, tt = (tti, 7r2, tts, 7r4) maps the P^s onto points, 

7ri(pi) = x, 7r2(pi)=7r3(pi) = 7r4(P^) = 0, (3.27) 

and is a birational isomorphism outside the P^s whose inverse is given by 

wi = TTi , w;2 = -T ( TTS + vr4) , Z = ; (3.28) 

^ TTs + 7r4 

By taking the difference between (3.25) and (3.26), we obtain the following 
algebraic constraint 

Mo : 7r| = ttI + vri - ^^(7^^) (3.29) 
that defines the singular Calabi-Yau geometry. 

Example 2: The resolved geometry for the model (2.18) takes the form 

z' = 1/z , w'i = z~^wi , w'2 = z^W2 + zQ{wi/z) + z'^P{w\) — w\z . (3.30) 

Prom (3.7) and (3.10) we see that the P^s sit at 

wi{z)=x + yz, w[{z') = y + xz' , (3.31) 
Pix)-Pix + yzl^ Qiy + xz')-Qiy)_ 

z z' 

The blow down map can be constructed by trial and error. For example, 
starting with the ansatz 

TTi = —z'w2 + • • • = —z'^W2 — Q{wi/z) + Wi — zP{wi) + • • • , (3.33) 

we see that the missing part can be adjusted to cancel the non-holomorphic 
piece Q{w\/ z) in the right hand side, finally yielding 

TTi = -Z'W2 + Q{w[) =Wl- zP{wi) - Z^W2 ■ (3.34) 



638 



PLANAR DIAGRAMS AND CALABI-YAU SPACES 



One can construct similarly 

1^2 = ZW2 + P{wi) = w'l — z'Q{w'i) + z''^W2 . (3.35) 

The construction of ir^ and 7r4 requires further thought. A hint is given by 
calculating tti and 7r2 on the spheres (3.31), (3.32), 



'^)=x, 7r2(P^)=y. (3.36) 



Consistently with the fact that the P^s are mapped onto points, the result is 
independent of z or z' . The result also suggests that -P(vri) and Q{tt2) could 
be useful globally holomorphic objects to consider. One then finds that 

TTS = -Z'^W'2 + z'''Qiw[) + z'PiTTl) - Z'w[ = -^(^l) " _ ^ (3.37) 

= z'w2 + z^P{wi) + zQ{7r2) - zw, = ^^""'^ '^^"^'^^ + w'2 , (3.38) 
have all the required properties. In particular 

7r3(pi) =7r4(P^) = (3.39) 
are z-independent, and the inverse is given by 

7r27r4 (7r2 - P(u;i))7r3 774 , . 

Wl=TTl + —— , W2 = > z = —— , (3.40) 

Q[-K2)-T^i P{tti)-tt2 Q(vr2)-vri 

with similar formulas for the primed variables. The singular Calabi-Yau is 
then found to be 

A^o : 7r37r4 = (P(7ri) - 7r2) {Q{'K2) - tti) . (3.41) 

Example 3: We consider the superpotential 

W{X, Y) = XY^ + V{X) + YU{Y^) . (3.42) 

Note that a term in could be generated by a simple shift in X. It is 
useful to introduce 

V'{x) = Fix) = -F,{x^) - xF2{x^) , 
Q{y) = U{y^) + 2y^U\y^) = -G{y^). ^ ' ' 

The resolved geometry is 

z' = 1/z , w[ = z~'^wi , W2 = z^W2 + zQ{wi/z) + z^P{wi) + wl , (3.44) 
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and the equations for the spheres are 



wi{z) = x + yz, w[{z') = y + xz' , 



(3.45) 



P{x) - P{x + yz) , , 2 , Q{y + xz') - Q{y) 
W2{z) = ) W2[z ) = X -\ — (3.46) 

The blow down map is constructed using the same tricks as for Example 1. 
A special case of the map also appears in [29] and [12]. We find 



TTl = wj- zG{tT2) + z'^P{wi) + Z'^W2 = W2 + 



G{w'i) - G{7T2) 



vr2 



7r4 



= —ZW2 — Piwi) = w'1 + z'Q{w'i) — z'^W2 , 

= 2;'7r2 - w'yF2{t\\) - z'F\{'k{) 

, F^{wl) - FiiTTi) + wi{F2{wl) - F2(7ri)) 

= —W2 H 5 

z 

= w'iir2 — -2;'7riF2(7ri) — w'iFi{Tri) 

, wi{Fi{wj) - Fi(7ri)) + wlF2{wj) - 7riF2(7ri)) 
= —W1W2 H 



The blow down can be inverted using, for example, the relations 
, (7r2 - Fi(7ri))7r3 + F2(7ri)7r4 



w'2 



(7r2-Fi(^i))2-7riF|(^i) 
7ri7r3F2(7ri) + (1x2 - Fx{Ti\))TXi 
(7r2-Fi(7ri))2-7riF|(7ri) 
, G(7r2) - Giw'l) 

TTl H ; 



(3.47) 
(3.48) 

(3.49) 
(3.50) 

(3.51) 
(3.52) 

(3.53) 



The P^s are mapped onto points, 
7ri(P^) = x^ , 7r2(P^) = , 7r3(P^) = -yF2{x^) , 7r4(P^) = xyF2{x\ (3.54) 
and the singular geometry is given by 



: 7r2 



(7r2-Fi(7ri))'-7riF|(7ri)" = 

7r| - TTiTrl - G(7r2) \{t^i - Fi(7ri)) tts + 7r4F2(7ri)l . (3.55) 



3.3 The singular geometry and privileged coordinates 



The singular points in a geometry 

A^o : ^o(7ri,7r2,7r3,7r4) = 



(3.56) 
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are obtained by solving the equations 

ao = , d^o = (3.57) 

simultaneously. By construction, the points 7ri(P^), which are the images of 
the two-spheres, are singular. They are associated with the one dimensional 
representations of the algebra of classical equations of motion, since the 
P^s sit at the extrema of the superpotential as a function of commuting 
variables. A deep consistency requirement is that there are also singular 
points associated with the higher dimensional representations [12], because 
the geometry encodes the most general solution of the matrix model. The 
full set of representations must thus be found from the algebraic equations 
(3.57) over commuting variables. This is possible if there is a correspondence 
between some combinations of the coordinates and the Casimir operators 
that characterize the representations. These combinations are related to the 
privileged coordinates that we have already discussed in Section 2.3. There 
can also be additional singularities, that are reminiscent of the commonly 
found branch cuts on unphysical sheets for the resolvents. 

Example 1: The singular points of the geometry (3.29) are at 

P(7ri) = , TTs = TTs = 774 = , (3.58) 

which yields the obvious identification 

X = TTi = X . (3.59) 

Example 2: For (3.41), equations (3.57) are equivalent to 

TTi = Q{t^2) , 7r2 = f (vTl) , TTs = 774 = . (3.60) 

This is perfectly consistent with (2.20), with the identification 

X = TTl=X, Y = TT2=y, (3.61) 

that is also suggested by (3.36). 

Example 3: The algebra of classical equations of motion associated with 
the model (3.42), 

A: {X,Y} = G{Y'^), y2 ^Fi(X2)+XF2(X2), (3.62) 

is very similar to the one studied in Section 2.2.3, and provides an example 
with both one and two dimensional representations. The generic singular 
points of the associated M.o (3.55) are found to satisfy either 

4771772 = G2(772) , (^2 " Fi{lTl)f = 77iF|(77i) , 

(Fi(77i) - 772)G(772) ^r(^\V(^\ Ci R-^^ 

773 = — ) 774 = -G(772ji'2(77lj , (3.63) 
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which yield the one dimensional representations, or 

i^2(7ri) = 0, 7r2 = Fi(7ri), tts = 7r4 = , (3.64) 

which yield the two dimensional representations. We obtain the identifica- 
tions 

= TTi = x^ y2 ^ 7r2 = y^ (3.65) 
which also make (3.63) and (3.64) perfectly consistent with (3.54). 



3.4 The deformed geometry 

Singularities of A^o lie at the classical critical points of the superpotential. 
When the coupling S is turned on, the eigenvalues spread and the critical 
points are replaced by cuts. The defining equation (3.56) of ^Ao is then 
deformed in such a way that the singularities are replaced by three-spheres. 
The equation of the deformed space is of the form 

M : ^(7ri,7r2,7r3,7r4) = 

^o(7ri, 7r2, TTs, 1:4) + S ^ Cabcd KA'^l'^i = • (3-66) 
(a,6,c,(i)eN4 

To find the allowed monomials vr^TrgVrgTr^, we must impose that the only di- 
vergences that occur in the period integrals (2.34) are either 5-independent 
or linear in 5' and logarithmic. This constraint comes from the renormaliz- 
ability of the gauge theory on the brane. A logarithmic divergence propor- 
tional to S is absorbed by the standard renormalization of the gauge coupling 
constant. It is important to realize that divergences proportional to higher 
powers of S are not allowed, even if they are logarithmic. For example, this 
subtlety must be taken into account to find the correct deformation of (3.55). 
Another constraint, from analyticity, is that a monomial that is generically 
forbidden cannot appear in a special case, for example when a coupling in 
the matrix model potential is set to zero. After all those constraints have 
been taken into account, and up to coordinate redefinitions, the number of 
deformation parameters Cabcd must match the number of irreducible repre- 
sentations, because there should be a uniquely defined deformed space for 
each vacuum (2.3). 

It is usually straightforward to guess the general form of the deformations 
up to coordinate transformations. For example, the deformed version of 
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(3.29), (3.41) and (3.55) are respectively 

M: 7r| = 7r| + 7ri-p2(7ri) + 5A(7ri), (3.67) 
M 

M 



7r37r4 = (P(7ri) - tts) (Q(7r2) - tti) + 5A(7ri, tts) , (3.68) 



T^2 



(7r2-Fi(7ri))'-7riF|(7ri) +5A(7ri,7r2 



7r| - G(7r2)F2(7ri)7r4 - ttittI - [^(712) (7r2 - Fi(7ri)) + 5A(7r2)]7r3 , (3.69) 

where A and A are polynomials. The constraints on the degrees of the 
polynomials can be found by explicit calculations. For the one-matrix model 

(3.67) it is well-known that deg A = degP — 1. For the model (3.68), wc will 
find in the next Section that the most general monomial in A is 7r"7r2 with 
< a < deg P — 1 = — 1 and < 6 < deg Q — 1 = c?y — 1. The number of 
deformation parameters is thus dxdy, consistently with (2.20). In the case 
of (3.69), three cases must be distinguished. If do = degG = 0, A must be 
a constant and the allowed monomials in A are vr" for < a < d — 1 and 
7r"7r2 for < a < [{d — l)/2]. If = = 1, A is a constant and A is linear 
in 7r2 and independent of tti. In all the other cases, deg A = — 1 and 
A is a linear combination of terms 7r^7r2 with < a < d — 1 for < 6 < 
da ~ 1, < a < d - 2 ioT do < b < 2da - I, and < a < [{d - l)/2] - 1 
for b = 2dG- It is easily to check that the total number of parameters is 
given by max(d -|- 2, 2ddG) + [{d — l)/2], matching the number of irreducible 
representations of the algebra (3.62). 



4 Solutions from the geometry 



The aim of this Section is twofold. First we want to explain how to get 
geometrically the resolvents for various matrices. Second we provide full 
calculations for the models (2.18) and (3.42). The solution of this latter 
example was not known, but we will be able to check the results in a par- 
ticular case in the next Section. Let us emphasize that any model of the 
form (1.1), at least for two matrices, should be in principle solvable using 
the same strategy. 



4.1 The resolvents from the geometry 

The idea is that, since any multi-matrix model can be formulated as a one- 
matrix model by integrating out all but one matrix, the Calabi-Yau geometry 
for a general model should have some common features with the geometry 



FRANK FERRARI 



643 




(z,w) 



S2CS3 



x-plane 



b 




Vi=V2 



Figure 4: Integrating over S"^. We first integrate over 5^ C (left inset) and 
then over a meridian of (right inset) . Prom the matrix model perspective, 
the meridian [a, b] is either an eigenvalue-filled interval corresponding to a 
branch cut in Figure 1, or to a intervals of the type [bx, clr'] or ]/xo, bx]- 

(3.67) for the one-matrix model. The relevant property of the latter geometry 
is that it is a fibration of the deformation of the simplest ALE space C^/Z2, 



The space = + + \ has a single S"^ of holomorphic volume A. 
When \{z) = 0, the two-sphere shrinks. Classically this is equivalent to the 
equations of motion. The base coordinate z describes the fluctuations of S"^ 
and is associated with the matrix X, z = x as in (3.59). We thus expect 
that the more general spaces we have to deal with are natural fibrations 
over bases parametrized by the privileged coordinates discussed in Section 
3.3. The fiber over a point x is in general much more complicated than 
a simple deformed 'C? jTL^ space. For example, the multi-valuedness of the 
effective potential implies that the fibers must contain several iS^s. 

To make those ideas quantitative, we propose that integrating over the 
two-spheres in the fibers yields the discontinuity of the resolvent across a 
branch cut. Let us introduce the coordinate x associated with the matrix 
X, and let us denote by \x the interior product associated with the vector 
field d/dx. If is the resolvent for X on the physical sheet, and g^ is the 
analytic continuation of g^ through a branch cut, then 



This is the fundamental equation that unambiguously determines g , as can 
be seen for example by taking x to be in the support of and using (2.10) 



(4.1) 




(4.2) 



and (2.11). 
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We can now integrate over a three-sphere in the full three-fold geometry 
by first integrating over the S"^ in the fiber and then over a meridian [a, h] 
of in the base, as depicted in Figure 4. Equation (4.2) implies that 

(f n= [' [ n= [\g''{x)-g''ix))dx= ({g''{x)dx, (4.3) 

JS^ Ja JS^CF^ Ja J 

which links the special geometry relations in the Calabi-Yau (2.34) and ma- 
trix model (2.12), (2.17) forms. 

As we have said, there are in general several S'^ C F-x over which to 
integrate. The relevant S'^s can be easily identified by looking at the classical 
5^0 limit of (4.2). Equations (4.2) and (2.14) imply 



lim / ix^ = fh,ci{x) , (4.4) 

where the background force was defined in (2.16) (for the D > 2 dimensional 
representations, one must take into account the terms discussed in Section 
2.2.3 in the classical background force). Note that a rigorous consistency 
condition that follows from the second equation in (2.34) is that the non- 
compact integrals of lim^^o fs'^cF. related to the critical values 
of the potential W. We believe that the stronger condition (4.4) holds if 
and only if the ansatz (4.2) is correct. Equation (4.4) is very handy to 
fix the overall normalization of J7. Let us also point out that the classical 
equations of motion, including the higher dimensional representations, are 
strictly equivalent to the condition ci(^) = 0) 



trdW" = 0^1im/ ix^ = . (4.5) 

The basic requirement of renormalizability is equivalent to the fact that 
the quantum corrections to jg2f-p ^x^ must vanish faster than 1/x when 
X — > 00, except for terms linear in S that may go as 1/x. As explained in 
Section 3.4, this must restrict the number of moduli of the deformed space 
to be equal to the number of irreducible representations of A. 

EXAMPLE: Let us consider a non-compact Calabi-Yau threefold in defined 
by the equation 

M : = f{z,w)v'^ + g{z,w)v + h{z,w) , (4.6) 

where /, g and h are polynomials. The nowhere vanishing holomorphic 
three-form is 

^^^ d.AdzAd^ ^ (4.7) 
2i7ru ^ ' 
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where is a normalization constant. We assume that the privileged coor- 
dinate X is expressed in terms of the variables z and w only, so that 

dv 

ix^ = -N- Aix(dzAdw). (4.8) 

ZITTU 

As in Figure 4, the two-spheres arc generated by closed contours encircling 
the branch cut [i^i,i'2] in the u-plane. Equation (4.2) then yields 



\/{v -vi){v -V2) \/f{z, w) 

._^rw^. (4..) 

The North and South poles N and S lie on the algebraic curve vi = V2, 

C : g^{z,w) -4f{z,w)h{z,w) =0. (4.10) 



4.2 The two-matrix model with XY interaction 



By introducing the privileged coordinates (3.61) z = ni = x and w = 1^2 = y 
associated with the matrices X and Y respectively, the geometry (3.68) can 
be put in the form 

M: u^ = v'^ + C{x, y)=v^ + {P{x) - y) {Q{y) - x) + SA{x, y) . (4.11) 

Equation (4.9) yields 

5^(x)-5^(x) = -iV / dy. (4.12) 

Let us denote by y = yi{x) the + 1 solutions to the equation 

C(x,y(x))=0, (4.13) 

with labels chosen in such a way that, classically, yi,c\{x) = P{x) and 
yi,c\{x) = Q~^{x) for i>2. Consistency with (4.4) and (2.21) then requires 
that = ±1 and 

g''{x)-g''{x)=yi{x)-y^{x). (4.14) 

Let us note that the physical sheet branch cuts of the resolvent correspond 
to and only to yi = y^, z > 2. Prom (4.14) we thus deduce that and — yi 
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have exactly the same branch cuts and discontinuity across branch cuts, and 
thus must be equal up to some entire function. From (2.9) we finally get 

9^(x) =P(x)-yi(x). (4.15) 

The same analysis could be repeated to compute the resolvent for the matrix 
Y. 

Let us now implement the constraint of renormalizability, which restricts 
the large x behaviour of yi{x) — yi{x). These restrictions actually apply to 
the S-dependent part of yi{x) and yi{x) separately, because generically no 
cancellation can occur. A straightforward calculation then shows that the 
allowed monomials x^y^ in A(x,y) are such that 

a + hdx <dxdY -I, a + h/d^ < dx + {dy - l)/dY - I ■ (4.16) 

This is equivalent to a < dx — 1 and 6 < (iy — 1, showing that there are 
exactly dxdy free parameters in A, matching the number of irreducible 
representations of A. 

Summary: If C{x, y) = {P{x) - y) {Q{y) - x) + SA{x, y), with 

dx — 1 dY ~ 1 

A{x,y)= J2 E (4.17) 

then the resolvents g-^ and for the model (2.18) satisfy 

C{x,P{x)-g''{x))=C{Qiy)-g^{y),y)=0, (4.18) 

with the asymptotics (2.9). The dxdy parameters Cij are determined by a 
choice of vacuum through the dxdy conditions (2.12). 

The above result is in perfect agreement with the known solution from 
the loop equations [14], which yields 

^ \n x-X y-Y I ^ ' 

Note that in the loop equation approach the vacua are naturally character- 
ized by a set of basic correlators, whereas in the geometric approach the 
natural parameters are the filling fractions (2.5). There is a one-to-one cor- 
respondence between those two sets of parameters. 
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4.3 The generalized Laufer's matrix model 

We now solve the model (3.42), with deg V' = d and deg U = da- When U = 
and V{x) = x'^, the relevant Calabi-Yau (3.69) reduces to the geometry 
originally constructed by Laufer in [29]. With an obvious redefinition of 
coordinates, (3.69) can be cast in the form (4.6), 



M : = zv^ + 



{w - Fi{z)) G{w) + Sk{w) 



+ -G^{w)F^{z)+w\{w-Fi{z)f -zF^{z)\ +SA{z,w), (4.20) 
and the algebraic curve (4.10) is given by the equation 

C{zM = {zw - ^) {{w - F,iz))' - zFiiz)) 

+ SzA{z, w) -^S{w- Fi{z)) G{w)A{w) - ^S^A'^{w) = . (4.21) 

4.3.1 The resolvent for X 

By using the privileged coordinates (3.65) z = vri = and (4.9) we get 

g^{x) - g^{x) = -2NAw{x) , (4.22) 
where Aw{x) denotes the difference between two roots of the equation 

C{x'^,w{x)) =Q. (4.23) 
There are max(3, 2 + 2c?g) roots, labeled such that classically 

^i,ci(a^) = -y'{x) , W2,c\{x) = -V'{-x) , 

Ax'^Wi,a{x) = G^Wi,a{x)) forz>3. (4.24) 

The classical background force for one dimensional representations is simply 
/b,ci(a^) = —y'^ix) — V'{x) with 4a;^?/^(x) = G^(y^(a;)). Consistency with 
(4.4) and (4.5) then immediately requires that N = ±1/2 and 

(x) - g^ (x) = wi (x) - Wi{x) , i>2, (4.25) 

where the cases i = 2 and i > 3 correspond to two and one dimensional rep- 
resentations respectively. Finally the branch cut structure and asymptotics 
of the resolvent imply 

g^{x)=wi{x) + V'{x). (4.26) 



648 



PLANAR DIAGRAMS AND CALABI-YAU SPACES 



Let us now study the normalizability constraints that follow from (4.25). 
One must be careful that the term in 5^ in (4.21) does not produce any di- 
vergences, while terms in S may yield logarithmic divergences. The relevant 
constraints on the monomials z^'nfi in /S.{z,w) read 

2a + bd < max(d + 2,2ddG) + 2[(d- l)/2] -2, (4.27) 
(2dG -l)a + b < {2dG - l)d - do for do > 1 and d > 2 if = 1 , (4.28) 
a + b < 1 for d = dc = 1 , (4.29) 

and the relevant constraints on A{w) read 

2ddegA < max(d + 2,2ddG) + 2[(d-l)/2] -1, (4.30) 
ddegA < max(d+2,2ddG) + 2[(d-l)/2] -d(l + dG). (4.31) 

The inequalities (4.27), (4.30) and (4.31) are obtained by studying the 
asymptotics of the roots wi or W2, and (4.28) and (4.29) follow from looking 
at Wi for i > 3. It is straightforward to check that the general solution to 
this set of inequalities yields exactly the deformed geometry described at the 
end of Section 3.4. 

Summary: Let C be defined by (4.21), with the constraints on A and A as 
in Section 3.4. The resolvent g-^ for the model (3.42) satisfies the degree 
max(3, 2 + 2dG) algebraic equation 

C{x\g''{x)-V'{x)) =0, (4.32) 

with the asymptotics (2.9). 

4.3.2 The resolvent for Y 

Using (3.65), (4.9), and the normalization constant N deduced in the pre- 
ceding subsection, we obtain 

/(y)-5''(y) = 2yAx(y), (4.33) 

where Ax{y) denotes the difference between two roots of the equation 

C{x{y)^y^)=0. (4.34) 

Let us separate the 2(d-|-l) roots into two sets {xi} and {— Xj} characterized 

by 

xi,ci{y) = ' V [xi^M) + = for ^ > 2 . (4.35) 
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One and two dimensional representations correspond respectively to xi = 
—Xi and Xj = —xi for i,j > 2. From the classical background force in one 
dimensional representations, /b,ci = G(y^) — 2yx{y) with + V'{x{y)) = 0, 
we deduce that 

/(y) - f{y) = -2y {xi{y) + xj{y)) , (4.36) 

where i > 2 and j > 1. The resolvent is the only analytic function 
satisfying (4.36), with no other branch cuts than those associated with the 
representations of A, and that goes as S/y at large y on the physical sheet. A 
naive guess might have been to identify with one of the root Xj for z > 2, 
but this cannot work since for example there would be unphysical branch 
cuts corresponding to the permutation of the root Xi with another root xj. 
It is easy to eliminate those latter branch cuts by taking a permutation 
invariant sum, normalized in such a way that (4.36) is satisfied, 

d+l 

9^{y) = -y^Xi{y) + polynomial . (4.37) 

i=l 

The only possible remaining unwanted branch cuts in Xi would corre- 
spond to the permutation of xi with — .xi, but it is actually straighforward 
to show that there are no cuts permuting any root Xi with its opposite. In- 
deed, the associated 5 = double points do not open up because the curve 
(4.21) factorizes at z = for any S, 

C(0,u;) = -^ (w-Fi{{)))G{w) + SK{w) ^ (4.38) 

Finally, we fix the polynomial part in (4.37) by looking at the large y asymp- 
totics. 

Summary: Suppose that V'{x) = Yfk=o''^k+ix'^ ■ The resolvent for the 
model (3.42) is given by 

d+l 

a'iy) = -yJ2^iiy)--Giy')-^ iid>2, (4.39) 

= -y {xi{y) + X2iy)) - \G{y') - ^^^^ + if d = 1 ,(4.40) 

where the sum is taken over half of the roots of the equation (4.34) satisfying 
the conditions (4.35). 

The resolvent 5^, being a sum of algebraic functions, must itself be an 
algebraic function. It is not difficult to find the algebraic equation satisfied 
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by g . Let us introduce the symmetric polynomials (Jq = 1 and ak = 
J2^ii " ' ^ik ioT 1 < k < d + 1. The resolvent is essentially —yai. We can 
write 

d+l 

C{x\y^) = {-l)%+^y^\{{x-Xi{y)){x + Xi{y)) 

i=l 

= {-Iftl^.y' J2 i-'^r<^k^k'x"'^'-'-'' ■ (4.41) 

0<A;,fe'<d+l 

Comparing with (4.21), we get d+l quadratic equations for the d+l unknown 
ak- By looking at the constant term and the classical limit (4.35), wc find 
that aa+i = (-1)'^+^ {G{y^)iti + y^) + SA{y^)) /{2ytd+i). There remains 
d quadratic equations for ai,...,ad- Geometrically, this means that the 
resolvent lies at the intersection of d quadrics. By elimination of variables, 
we can find a degree 2'^ algebraic equation satisfied by . For example, 
when d = 1, we recover the hyper elliptic curve of the ordinary one-matrix 
model obtained by integrating out X. 



5 Solution from the loop equations 

The general class of models (1.1) ought to have some very special mathe- 
matical properties that make possible a solution in terms of an algebraic 
variety. In some sense, the loop equations should be integrable in the planar 
limit. A general analysis is far beyong the scope of the present paper, but 
we have been able to solve the example (2.23), which is a particular case of 
(3.42) for G = —a. Wc find a precise match with the result of the previous 
Section, providing a very non-trivial test of the geometric approach, and, 
most importantly, of the underlying conjectures on which it is based. The 
test is particularly stringent because the model (2.23) is associated with a 
genuine J\f = 1 gauge theory that is not a deformation of an underlying 
J\f = 2 theory. In particular, the solution cannot be found by perturbing a 
Seiberg-Witten curve, unlike all the examples studied so far in the literature. 

5.1 General loop equations 



We consider the matrix integral 



(5.1) 
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for an arbitrary polynomial 

d+l 

V{X)=Y,tkX^lk. (5.2) 

fc=0 

The basic object we are going to compute is the generating function 

T{x,y) = s(^--^^). (5.3) 

\n X — X y — Y / 

It can be expanded, either at large y or at large x, 

V{x,y) = Y,g^{x)y-'-' =Y,9k {y)^'"-' , (5-4) 

fe>0 A;>0 

in terms of generalized resolvents 

The ordinary resolvents are = and g^ = g^ ■ 

Let us now consider the following variations in (5.1), 

5X = , 5Y = — , (5.6) 
X — X 



5X = -(— ^— + — (5y = 0, (5.7) 

2\x-Xy-Y y-Y x-X) ' ' ^ ' 

SX = 0, 6Y = U^ ^ + ^^ ^ '—). (5.8) 

2\x-Xy-Y -x-X -x-Xy-Yx-XJ ^ ' 



The remarkable property of these changes of variables is that, due to many 
cancellations of terms, the associated loop equations close in the planar limit 
under the correlators of the form (trX^y). After calculating the jacobian 
of the transformations, taking into account the variation of the classical 
potential, and using the factorization of multi-trace correlators at large AT, 
we get 

(5^(x) -y^- V\x)) r(x,2/) = -(y - ^)5^(x) - 5 

'tr 1 V'{x)-V\X)\ 



V 1^ /' (^-9) 

n y — Y X — X I 

V{x,y)V^-x,y) = ^^(x) + ^^(-x) - y (r(x, y) + r(-a;, y)) 

-^(r(^,2/)-r(-a^,2/)) • (5-10) 
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5.2 The resolvent for X 



The functions , and in particular the resolvent for X, are obtained by 
expanding the loop equations at large y. Prom (5.9) we get the following 
recursion relations, 



5f(.) = si^i^-f ,^(x), (5.11) 
nx 2x 

9k+2{x) = {g''{x)-V'{x))g^{x) + SAk{x), (5.12) 

where the Ak{x) arc degree d — 1 polynomials that can be computed from 
g^ and that are defined by 

/tr V'(x)-V'iX) \ 
A,(x) = (-y ). (5.13) 

The relations (5.11) and (5.12) determine all the g^s, and thus r(x,y), from 
the single function g'^{x). Moreover, from (5.10) we get 

9q+2ix) + 4+2{-x) = {gf+lix) - gf+ii-x)) - Yl 9kix)9k' i-x) ■ 

k+k'=q 

(5.14) 

By using (5.11) and (5.12) in (5.14) we obtain, for each q, equations that 
close under g^ (x) and g^{—x). Any two such independent equations then 
determine unambiguously g^{x). It turns out that the equation for g = 1 
does not yield anything new, so wc use g = and q = 2. Eliminating g-^{—x) 
we get, after a lengthy but straightforward calculation, a closed equation for 
g^{x) only. It is a cubic algebraic equation of the form 

C{x^,g^{x) - V'{x)) = 0, (5.15) 

where C is defined exactly as in (4.21) for G = —a. The deformation poly- 
nomials are found explicitly in terms of correlators to be 

A{w) = 2(^y), (5.16) 

A(x^ w) = A2(x) + A2(-x) + xF2{x'^) (Ao(x) - Ao(-x)) 

-Fi(x2) (Ao(x) + Ao(-x)) + ^ (Ai(x) - Ai(-x)) 

+ (Ao(x) + Ao(-x)) w . (5.17) 

The form of A and A is exactly as discussed at the end of Section 3.4: the 
term independent of u; in A is a polynomial of degree d — 1 in z = x^, 
the term linear in w is a polynomial of degree [{d — l)/2] in z, and A is 
a constant. We have thus found a precise match with the result obtained 
from the geometric approach. Note that similar cubics were found from loop 
equations in [20] for models that are deformations oi M = 2 gauge theories. 
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5.3 The resolvent for Y 



One way to compute the resolvent for Y is to expand the loop equations at 
large x. Prom (5.9) we get a set of linear equations 

E tk+i9L,{y) = SC^ Y)6,,o - ^ (^X^-i) (1 - 5,,o) + Sy {'^X") 

It f i ft 

k=0 

-y'gy{y) + S ^ {^^X>^) g^iy) , (5.18) 

k+k'=q-l 

and from (5.10) we get a set of quadratic equations 

yd,-i{y) = s C^x^'-') - f 52^2(2/) + ^ E {-^)'9l{y)9l'{y) ■ 

k+k'=2q-2 

(5.19) 

The 2d equations (5.18) for < q < d—1 and (5.19) for 1 < q < d close under 
the 2d unknown for < k < 2d — 1. We can for example express linearly 
the 5^ for cZ < A; < 2(i — 1 in terms of the for < A; < d — 1 by using 
(5.18), and then obtain a set of d quadratic equations for the d unknown 
g^ , . . . ,gj_i from (5.19). Wc thus discover that the resolvent for Y lies at 
the intersection of d quadrics, consistently with the result of Section 4.3.2. 
However, the direct calculation showing that the set of quadratic equations 
obtained from (5.18) and (5.19) on the one hand and (4.41) and (4.21) on 
the other hand are equivalent turns out to be extremely tedious. 

We are thus going to provide a much simpler proof using the fact that we 
have already computed the resolvent for X in Section 5.2. Let us consider 
the equation for the unknown x{y) 

5^(x(y))-y'(x) = y2. (5.20) 

We know that x{y) then automatically satisfies (4.34). Not all the solutions 
to (4.34), though, satisfies (5.20). This is best seen by plugging (5.20) in 
(5.9), which yields 

a\/9 / ^,tr 



,,.,,)(5^nfMzp^)).„. (.21) 

\ n y — Y x[y) — X I 

This is a degree d + 1 algebraic equation for a;(y), and thus only half of the 
roots of (4.34) can actually satisfy (5.20). By taking the classical limit 5 — > 
of (5.21), we see that those d + 1 roots are precisely the roots Xi that we 
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have used in Section 4.3.2 and that were characterized by (4.35). By looking 
at the coefficient of x'^'^^ (for the overall normalization) and of x'^ in (5.21), 
we then immediately get the sum a\ = Yli=i ^iiv)^ 

a, = __^+ « if d>2 (5.22) 

U+i 2y y 

= + if d=l. (5.23) 

Those equations are equivalent to (4.39) and (4.40). Moreover, we see that 
the higher symmetric polynomials Ck, 2 < k < d, are related to the general- 
ized resolvents iov k < d — 1. Since the degree d is arbitrary, we conclude 
that the Calabi-Yau geometry actually encodes the full generating function 

r. 



6 Discussion 



The geometric approach to matrix models is singled out by its aesthetic fea- 
tures and its relationship with some of the deepest insights in gauge theory 
and string theory. It provides an entirely new perspective on the problem 
of summing planar diagrams, and suggests that a whole new class of mod- 
els could be solved. The examples that we have studied show that, when 
appropriate, the geometric approach, which reduces to the calculation of a 
blow down map, is much more powerful than standard techniques. There 
remains, however, many open problems, some of which we review below. 

Matrix model technology: it is now clear that all the irreducible rep- 
resentations of the algebra A of equations of motion should be taken into 
account, even though most of the classic matrix model literature considers 
only one dimensional representations, often with the additional one-cut re- 
quirement. We have demonstrated in Section 2 that new qualitative features 
occur for higher dimensional representations. It would be desirable to de- 
velop the analytic techniques to deal with the saddle point equations in those 
cases, and to have a more general derivation of the property of eigenvalue 
entanglement. Another basic problem is to understand the general condi- 
tions under which special geometry relations can be derived. The geometric 
approach suggests that Casimir operators play a special role in this respect. 
Another interesting aspect is that the difficulty of a given model seems to 
be directly related to the complexity of the algebra A (the representation 
theory and the structure of the center). It is not clear how this translates in 
the language of loop equations. An interesting concrete problem would be 
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understand what makes the algebra and the loop equations associated with 
the models (1.1) special. 

Algebraic geometry: the infrared slavery of gauge theories suggests that 
all the Calabi-Yau geometries of the form (3.16), or even (3.6) when the 
corank of the Hessian of the potential W at critical points is at most two, can 
be blown down. It is not clear how this works mathematically. In particular, 
it would be desirable to understand the relationship with the results of [28] 
on Gorenstein threefold singularities. A naiVe guess would have been that 
blow down can only be found for the cases constructed in [12], but clearly 
the models (1.1) are much more general, and already the model (3.42) that 
we have solved provides a counter-example. A startling fact is that the 
calculation of the blow down map vr is essentially equivalent to finding the 
solution of the associated matrix model. It might be possible to devise an 
algorithm that computes tt. Another interesting feature is the interplay 
between the non-commutative structure and the singularity structure of the 
blown down geometry, which must reproduce the representation theory of 
the model. In some sense the blow down map knows about D-branes. 

Geometric approach technology: we have explained in Section 4 how 
to calculate the resolvents from the geometry. It is not clear to what extent 
this approach can be applied in general. The issue is to understand the 
fibcrcd structure of the geometries, which is probably related to the structure 
of the center of A. We would also like to understand what arc the most 
general matrix model correlators encoded in the geometry, and also how to 
extract the non-planar contributions. The answer can probably be found in 
the topological string set-up [31]. 

D-BRANES ON Calabi-Yau: an interesting question is to ask what is the 
most general matrix model that can be engineered by putting branes in a 
Calabi-Yau. This is the "reverse geometric engineering" problem [32]. When 
there is no moduli space, which is the generic case we have been considering, 
we have very little insight into the solution of that problem. We have been 
able to construct in Section 3 a large class of models, but it might be possi- 
ble to find additional theories by considering more general perturbations to 
the geometry (3.2). For example, one would like to know, given a in a 
Calabi-Yau, what is the most general geometry for which the obstruction to 
the versal deformation space of the is integrable in terms of a potential, 
taking into account the non-commutativity of the variables. The solution is 
probably most naturally expressed in terms of some "regularity" conditions 
on the algebra A. Even if we could understand the reverse geometric en- 
gineering, there would remain the fundamental question of why and when 
the geometric transition conjecture is valid. The phenomenon it describes is 
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very reminiscent, both physically and mathematically, to the "continuation 
to negative radius" found in two dimensional a models (see for example [33]). 
In this latter case, we know that some non-geometrical phases are possible. 
This suggests that we may presently only catch a glimpse of the full story for 
what can happen beyond the singular (infinite gauge coupling) point. The 
matrix models provide in principle a powerful tool to study that question. A 
full understanding would amounts to describing the space of vacua in string 
theory. 

M = I GAUGE THEORIES: we have used throughout the full non-perturbative 
Dijkgraaf-Vafa conjecture. This includes, consistently with special geome- 
try, keeping the Veneziano-Yankielowicz term and the higher powers in the 
glueball superfields Si, even when the latter are perturbatively zero in the 
chiral ring. There is no proof of this conjecture at the moment. When a 
direct solution of the matrix model exists, as described in Section 5, our 
result can be interpreted as providing a non-trivial check of the consistency 
betweeen the Dijkgraaf-Vafa conjecture and the geometric transition picture. 
The two-matrix models we have studied are particularly interesting because 
they are not deformation oi Af = 2 theories (which requires either one or 
three adjoint fields), and thus there is no Seiberg-Witten curve. The matrix 
model is then the only tool at our disposal. It would be very interesting to 
work out the quantum space of parameters for these theories, along the lines 
of [9, 34]. In particular, it is in principle possible to study quantitatively 
the phase transition between the Higgs and confining phases on parameter 
space [35]. 

As a final comment, we would like to note that we have put the emphasis 
on matrix models and thus on theories with only adjoint fields. However, all 
the questions we have addressed are also relevant to the more general set-up 
of quiver theories. 

Acknowledgement s 

I would like to thank D. Berenstein, S. Katz, W. Lerche, C. Vafa and par- 
ticularly R. Dijkgraaf for very useful discussions, as well as the organizers of 
the Amsterdam summer workshop "String Theory and Quantum Gravity," 
June 2003, where some of the results of this article were first presented. This 
work was supported in part by the Swiss National Science Foundation. The 
author is on leave of absence from Centre National de la Recherche Scien- 
tifique, Laboratoire de Physique Theorique de I'Ecole Normale Superieure, 
Paris, Prance. 



FRANK FERRARI 



657 



A Saddle point equations and special geometry 

1 Continuum limit of singular sums 

Let us consider the sum 

n ^ X — Xi 
1=1 



We assume that when n oo, the distribution of eigenvalues (iCi)i<i<n goes 
to a smooth function p with compact support, and we want to compute 

a{x) = lim anix) . (A.2) 

n— >oo 

If X Support [p], we have 

a{x) = g{x) , (A.3) 
where the analytic function g{x) is defined by 

, , /•+~p(z)dz ^, 

If X € Support[/9], g{x) is ambiguous because g has a branch cut. Let us 
label the eigenvalues in such a way that xi < X2 < ■ • • < Xn, and pick 
Xj < X < Xj+i with 5+ = Xj+i — X and 5~ = x — Xj. Let us define 

?M = lim ^ • (A.5) 

6 [X) n-^oo Sn 

We can then compute 



(rx-S„ r+oo 
/ + / 

J-oo Jx+S+^ 



p{z) dz 



X — z 



1 (5+ fx) 

= -[g{x + i€)+g{x-ie)) +p{x)ln^^- (A.6) 

One can use the above formula to deduce equation (2.33) from (2.30). The 
same reasoning also yields the standard result 



n— »oo 



lim - — ; = - {g{xi + ie) + g{xi - ie)) (A.7) 



because at large n we have in that case = Xj+i — Xj = l/(np(xj)) 

Xi Xi—i S . 
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2 Analysis of the saddle point equations 



Let us discuss briefly the saddle point equation (2.33) supplemented with 
the condition (2.32). If the intervals are filled with eigenvalues in one 
dimensional representations, we have 



g^{x + ie)+g^{x-ie)+g^{-x) = V'{x) + ^ for x G . (A.8) 

Moreover, by taking the even and odd parts of (2.33) for two dimensional 
representations, and using (2.32) and (2.10), we get 



,2 



5- {x) 



Sp^'i-x) In Pf\ = ^ - Fi(a;2) - ^xF^ix") 



-\{9^{x + ie)+g^{x-ie)), (A.9) 



{x + ie) = {-X + ie) - 2xF2{x^) , for a; G U J'^ . (A.IO) 

Let us denote g^^^^ix) = g^ {x) on the physical sheet. The equation (A.8) 
shows that the cuts glue the physical sheet with a second sheet where 
the solution is 



Q,2 

4ond(a^) = - gLti-x) + V'{X) + ^ • (A.ll) 



Due to the term 5^g^(— x), this second sheet contains new cuts l'^ that arc 
the image with respect to a; = of the physical cuts Using again (A.8), 
we see that the new cuts glue the second sheet with a third sheet with 

5tWrd(^) = -5&st(^) 

-{-gLti-x) - gLtix) + v'{-x) + ^)+ y'{x) + ^ 

= 5&st(-^)-2xF2(x2). (A.12) 

Equation (A.IO) then shows that this third sheet is glued to the first sheet by 
the cuts and . Overall g^ has generically a three-sheeted structure, 
except when only two dimensional representations are present in which case 
the second sheet is absent. The general form of the cubic equation satisfied 
by g^ is derived both in Section 4 and in Section 5. A special form of this 
cubic in the one-cut assumption and for a = was found in [19]. Let us also 
note that models where a similar cubic appears have been studied recently 
in the literature [20]. 
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3 The effective potential for Y 

We have argued in Section 2.2 that the singularity structure found for the 
background force (2.29) was generic. We want to illustrate this fact by 
studying the effective potential for Y in the model (2.23), defined by 



I dXe-i^'-^^^'+^W). (A.13) 



We know that V^s must have singularities when deg F > 2 to make the 
associated saddle point equations consistent with two dimensional represen- 
tations. This is a little bit puzzling at first, because the effective potential 
for Y is essentially the same as the effective potential for the ordinary two- 
matrix model discussed in Section 2.2.2, except that it is a function of the 
square of the matrix instead of the matrix itself. A simple way to understand 
what is going on is to use the Itzykson-Zuber formula [5] to derive 



U.i<j {{vi - yj){yi + Vj)) 

i i<j 

This formula shows that the only singularities may be at Ui = yj or yi = —yj. 
To find which singularities actually occur, we can evaluate (A.14) in the limit 
S -^0. This yields 



^^^Ili<j{iyi-yj)iyi + yj)) 

where the function x(y^) = Y'~^ {—y'^) minimizes xy^+^(x)- When yi — >■ yj, 
the pole in the denominator of (A.15) is cancelled by a zero of the denomina- 
tor, consistently with the fact that no singularities are expected at y^ = yj. 
On the other hand, because x is d-sheeted, if we analytically continue yi from 
yi = yj to yi = —yj we can go to a different sheet xiyf) ~^ xiyf) ^^^d tl^us 
produce a pole in (A.15). We read from (A.14) that the relevant singular 
part in Ves{Y) is exactly the same as that for Vee{X) in (2.28). In partic- 
ular the eigenvalues will be entangled for two dimensional representations, 
implying that 

p^(y)=p^(-y) forye/^. (A.16) 



4 Derivation of special geometry relations 



An elegant way to derive special geometry relations is to set up a varia- 
tional formulation of the saddle point equations. Let us first deal with the 
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description in terms of the matrix X. One must be careful because 



dxi 

due to the logarithmic term in (A. 6). Let us consider the functional 
= -SV^sip^) + J dxdxp^{x)p^{x) In \x - x\ 

d+2 [{d-l)/2] 

+ Y,iK[ISp''{x)dx-SK)+ (KU_^Sp''{x)dx-SK) 
K=l "^^K K=l -^^K 



[(d-l)/2] 



+ fs{p^ix)-p^i-x))LKix)dx. (A.18) 

K=l •'^K 

By varying with respect to the Lagrange multipliers £k and £k we obtain 

/ Sp''dx=^(f g''dx = SK, f 5p^dx = -^/ g''dx = SK. 

(A.19) 

The notation for the contours is similar to that of Figure L The equa- 
tions (A.19) implement the relations (2.5), with Sk associated with two 
dimensional representations. By varying with respect to the Lagrange 
multipliers Lk{x), we get the constraint (2.32). Finally, by varying with 
respect to p'^{x) we get the saddle point equations 

= - f^^l^ + 2S fdx p^{x) In \x-x\+ for x G , (A.20) 
op yx) J 

= - J^^y + Lk{x) + 25 [dxp^{x)ln\x-x\+£K for x G I^, (A.21) 
op \X) J 

= -J^1^^Lk(-x) + 2S [dxp^(x)ln\x-x\ for x G . (A.22) 

op^ (x) J 

One can check that the derivatives of the above equations are equivalent to 
(A.8), (A.9) and (A. 10), with the identification 

L'k{x) = Sp^ {-x)ln^^ forxG/^. (A.23) 
S (x) 

Moreover, the partial derivatives with respect to Sk and Sk take very simple 
forms, 

lis-.-''" 
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Using (2.28) and (A. 20), it is straightforward to compute 

/ g''ix)dx = [ " i2g''' 



dx6p^{x). 



which yields explicitly 



= -^K + y^^-251n/xo + 0(l//xo), (A.25) 



lim (I g^dx-V(no) + Slnno]. (A.26) 



Similarly, using (A.21), (A.22) and (2.32) we get 



/ ff^ (x) dx = -Ik + //f , + . xT'"" . - 45 In /xo + 0(l//xo) , (A.27) 



'7k 

yielding 



= lim f/ 5^da;-y(/xo)-^^(-/xo) + 251n/xoV (A.28) 



Using the results of Section A. 3, we can repeat the discussion above 
by replacing the matrix X by the matrix Y. There is an additional K- 
independent term 

in (A. 24) because Ves{p^) is not linear in . This term does not affect the 
relations involving compact cycles, 

oSk OSk' J pi dSx dSK' Jpl^, 



B RG flow theorem for normal bundles 

In this Appendix, we provide a proof of the conjecture in Section 3.2.1 in 
the special case n = 1 relevant to two-matrix models. The proof is a direct 
consequence of two lemmas. 

Lemma 1: Let in the geometry (3.16) at a given critical point (x, y) of 

the supcrpotcntial (3.14). Then the transition function characterizing the 
normal bundle of the can be cast in the form 
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Proof: To calculate the transition function, we expand wi = wi{z) + 
W2 = W2{z) + 62, w[ = w[{z) + 5[ and ^2 = ^2(^) + ^'2 ™ (3-16), keeping 
only the linear terms in the Ss. By redefining 



S2^62-5,Y.z^''Y.^i-7'W/f-, 

«>3 j>0 



(B.2) 



we obtain 



i>0 j>0 



(B.3) 



Lemma 2: Let £^ be a holomorphic vector bundle over with structure 
group GL(2, C) and transition function 

T{z) = ( , f\ 2 . (B.4) 

^ ' \a + hz + cz'^ z"^) ^ ' 

Let r be the corank of the quadratic form Q = Then E = 0{r — 1) © 

0(-r- 1). 

Proof: We want to construct gauge transformations P^{z) and P'b,{z' = 
1/z) such that 

Wn^=("7 (B.5) 

The important property is that Pn{z), P-^^{z), Ps{z') and P^^{z') must be 
holomorphic functions. It is not difficult to check that the following matrices 
have the required properties: for ac — ^ and 6 7^ 0, 

n _fac/b — b — cz {a/b-z)z\ p _f -a? /b az'/b-l\ 

cz-b ) ' [b^ -ac + abz' {c-bz')z')' 

(B.6) 

for 6 = 0, a 7^ and c 7^ 0, 

p ^/-c(l + -z)/a l/c- z{z + l)/a\ 

(B.7) 



/2 



1 - az'jc z"/c - {z' + 1)1 a 
^ —a z' 

for a = a^,b = af3, c = a 7^ and /3 7^ 0, 

\-p^/a l-jjz/aj \ —a^ z -/3/a ' 
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for a = 6 = and cj^O, 



(B.9) 



for 6 = c = and a 0, 



= • -1/a) ' = (J _ , , 



(B.IO) 



and finally for a = 5 = c = 0, Pn = -Ps = 
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